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Abstract This paper reports my recent study"! 


on the shear viscosity of neutron-rich nuclear matter from a 


relaxation time approach. An isospin- and momentum-dependent interaction is used in the study. Dependence of 


density, temperature, and isospin asymmetry of nuclear matter on its shear viscosity have been discussed. Similar to 


the symmetry energy, the symmetry shear viscosity is defined and its density and temperature dependence are studied. 
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1 Introduction 


One of the major problems in nuclear physics is to 
understand the properties of nuclear matter under 
extreme conditions. This is related to the basic 
knowledge of the in-medium nucleon-nucleon (NN) 
interaction which in the present stage can still hardly 
be obtained from the ab initio theory of the strong 
interaction, i.e., Quantum chromodynamics. Our 
knowledge on the in-medium NN interaction today is 
mainly developed along two lines. In the first line, one 
starts from the bare NN interaction, which has been 
fitted very well from NN scattering data, together with 
phenomenological three-body interactions, so that the 
in-medium NN interaction and the properties of 
nuclear matter can be obtained through many-body 
theories. In the second line, the starting point is an 
effective in-medium NN interaction or Lagrangian, 
with the parameters fitted to the empirical nuclear 
matter properties obtained usually through mean-field 
approximations. 

Ten years ago, an isospin- and 
momentum-dependent mean-field potential (hereafter 
"MDI')_ was 
(especially the isospin effects) in intermediate-energy 


constructed to study the dynamics 


heavy-ion collisions together with an 
isospin-dependent Boltzmann-Uehling-Uhlenbeck 


(IBUU) transport model”. In addition to the good 


description of the empirical nuclear equation of states, 
the momentum dependence of this mean-field 
potential reproduces pretty good the optical potential 
extracted by Hama S et al. from elastic proton 
scattering data. The studies using this interaction 
have constrained the nuclear symmetry energy at both 


46] in 


subsaturation and suprasaturation densities 
addition to the dynamics of heavy-ion collisions, the 
MDI model has also been used in the study of 
thermodynamical properties of nuclear matter’”*!, It 
found that the 


momentum-dependent potential can be derived from 


was recently isospin- and 
an effective interaction with a density-dependent 
two-body interaction and a Yukawa-type finite-range 
interaction using Hartree-Fock calculation’, The MDI 
model thus serves as a useful effective in-medium 
interaction. 

In the past few years, the shear viscosity of 
the quark-gluon plasma (QGP) formed in relativistic 
heavy-ion collisions has attracted special attentions. 
From the study using a viscous hydrodynamical 


1°! it was found that the strong-interacting QGP 


mode 
behaves like a nearly ideal fluid, i.e., its specific shear 


viscosity is only a little 
[11] 


larger than the 
Kovtun-Son-Starinets boundary . Up to now large 
efforts have been devoted to study the shear viscosity 
of QGP" and hadron resonance gas!®”’! formed in 


relativistic heavy-ion collisions, while there are only a 
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few studies on the shear viscosity of nuclear matter 


formed in 
[20-24] 


intermediate-energy heavy-ion 


collisions . Even few studies are related to the 
isospin effects on the shear viscosity of nuclear 
matter'”*!, In this paper I will discuss my recent study!!! 
on the shear viscosity of nuclear matter using the MDI 
model mentioned above from a relaxation time 
approach, which gives an intuitive picture how the 
shear viscosity changes with the density, temperature, 


and isospin asymmetry of nuclear matter. 


2 Shear viscosity from a relaxation time 
approach 


The system concerned here is an isospin asymmetric 
nuclear matter with uniform neutron and proton 
density p, and p,, respectively, and the nucleons 
are thermalized with temperature T. The flow field 
u is static in the Z direction and its magnitude is 
linear in the coordinate x , ie. u,=cx and 
u, =u, =O. In the rest frame nucleons move with the 
flow field and follow Fermi-Dirac distribution n° in 
the equilibrium state. In the lab frame the equilibrium 
distribution is a simple boost by the flow field 
compared with that in the rest frame, denoted as n’. 
Due to NN collisions, the real distribution may be 
slightly away from the equilibrium distribution and is 
denoted as n, and the deviation from the equilibrium 
distribution ôn =n’ -n is much smaller than n°. 

The shear force between flow layers per unit 


area by definition can be written as 


TE (P. mu) pv) (1) 
In the above, z=n orp denotes the isospin degree of 
freedom, p,v, is the number of nucleons moving 
between layers per unit time per unit area, and 
p,—mu, is the momentum transfer per nucleon in the 
z direction. The nucleon velocity in the 
x direction v, can be further written as v, = p,/m} , 
with m, being the effective mass. Using the momentum 
distribution n, =n?+6n,to calculate the average and 
taking into account that the equilibrium momentum 
distribution n? is even in p,, Eq.(1) can be further 


written as 


F P d°p 
—= > d| (p, —mu,)—-6n, ; 2 
a7 de m "(2x7 K 


where d=2 is the spin degeneracy. 


In the following I will calculate ôn, by 
linearizing the isospin-dependent BUU equation as 
follows: 

ôn 5 

PAPO V, n (p) -V,U V ,n;(pi) 
fa a d'p, d’p, d'p, 
2)" (2x) (2x) (27) 

x[n, (p )n, (p,)A—n, (P-n, (p,)) 

=n, (pn, (p) -n, (P) -n,(p2))] 

x (27) 5° (B, + B- Pi- B;) ©) 

3 K OR À Jat 

aji Ee d al d P 
(2x) (22) (27) 
x[n, (p )n (p-n, (pD -n(p})) 
=n, (pn (P)0-n, (p )0-n (p) 


x(2z)° 5° (B, +P, -Pi — Py). 


à 


T,T 


2 


T,-T 


In the above, U,is the mean-field potential from the 
MDI modell, T is the transition matrix, the 
degeneracy d-—-1/2 takes the double counting of 
identical nucleon collisions into consideration, and 
1-n is from the Pauli blocking effect. Replacing n 


with n° in the first-order approximation, the 
left-hand side can be expressed as 

an,(p,) 

Soh yy "V,„n,(p,) -V,U, "Vpn, (Pi) 


| ôu, P,P, dn? ) a 
OX P dp P=Pi 
by using the properties of n°. Keeping only the 


on,(p,) term, the right-hand side of Eq.(3) can be 
expressed as 6n,(p,)/t,(p,), where r,(p,) is the 


relaxation time, i.e., the average time between two 


collisions for a nucleon with isospin 7 and 


momentum p; , and it can be written as 


: = = + = ? (5) 
T(P) TP" (p) T; (p) 


where r¥™*®(p )is the average time for a nucleon 
with momentum p, to collide with other nucleons of 
same (different) isospin, and they can be calculated 
respectively from 
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1\(2") g a 
Sine =|d p-dp,d cos6,,d cos 0 
(p,) í 2 Gn 2 2 12 


Ti 
x do, Pr Pa (6) 
dQ |m, (p) m,(p,) 
x[n; (p,)—n?(p,)n? (pi) 
=n; (pa)n; (p3) +n? (pin; (p51, 
1 (2x7 
- =d p;dp,d cos 6,,d cos 0 
7" (p) (2r) f 24P2 12 
do, P, PA (7) 
dQ |m, (p) m_,(p,) 


x[n? (p,)- n°, (p,)n? (pi) 
-n° (p )n?,(p}) +n? (pi)n?,(p})]. 


In the above 6,, is the angel between p, and p,, 
and @ is the scattering angel between the total 
momentum and the relative momentum of the final 
state. In free space the pp and np scattering cross 
sections are isotropic and they can be respectively 


parameterized as”! 


o =13.73—15.04/v +8.76 /v° +68.67v*, (8) 


pp(nn) 


On = —70.67 -18.18 / v + 25.26/v" +113.85v, (9) 


where the cross sections are in mb and v is the 
velocity of the projectile nucleon with respect to the 
fixed target nucleon. This parametrization describes 
very well the experimental data for the beam energy 
from 10 MeV to 1 GeV"!, It is worth to note that in 
the most probable collision energies the np scattering 
cross section is about three times the pp scattering 
cross section. In nuclear matter, the in-medium NN 
sections 


scattering cross are modified by the 


in-medium effective mass in the form of”! 


a 32 
medium U 
oan = Onn = ) , (10) 


where yy ( yy ) is the free-space (in-medium) 
reduced mass of colliding nucleons. 
Once the relaxation time 7,(p) is known, 
ôn, (p) can be calculated from 
0 
n, (p) =r, (p) Z PP ee 
ox p dp 
Using the definition F/A=-nņ(ôu,/ðx), the shear 
viscosity can be calculated from Eqs.(2) and (11) in 


(11) 


terms of the local momentum distribution n* as 


2.2 he 3 
p;p. dn, d°p 
n= 27d] DA (12) 
2 J pm, dp (27) 
by setting the magnitude of the velocity field to be 
infinitely small. Note that from Eq.(12) the shear 
viscosity is related to the local momentum distribution 
near the Fermi surface. 


3 Results and discussion 


Figure 1 displays the density, temperature, and isospin 
dependence of the relaxation time. In neutron-rich 


same 
n 


nuclear matter, r“ is larger while r®™° is smaller 
compared to that in symmetric nuclear matter as a 
result of less frequent np collisions and more frequent 
nn collisions. For the similar reason, 7;"" 
* is smaller compared to that in symmetric 


nuclear matter. From Eq.(5), the total relaxation time 
diff 


is larger 
while r 


is determined by r™ which is always smaller than 


t™! due to the larger np cross section than pp(nn) 
cross section in the most probably collision energies. 
Thus, neutrons have a larger relaxation time than 
protons in neutron-rich nuclear matter. It is seen in 
Panel (d) that the relaxation time decreases with 
increasing temperature due to more frequent collisions 
at higher temperatures. In addition, at lower 
temperatures the relaxation time peaks around the 
Fermi momentum, indicating a strong Pauli blocking 


effect for nucleons near the Fermi surface. 
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Fig.1 (Color online) Panel (a), (b), (c): Relaxation time for 
neutrons and protons as a function of nucleon momentum in 
symmetric (6 = 0 ) and asymmetric (5 = 0.5) nuclear matter at 
saturation density and temperature T =50 MeV; Panel (d): 
Relaxation time as a function of nucleon momentum in 
symmetric nuclear matter at different densities and 
temperatures. 
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Fig.2 (Color online) Density and temperature dependence of 
the shear viscosity ((a), (b)) and specific shear viscosity ((c), (d)) 
for symmetric (6=0) and asymmetric (6=0.5) nuclear 
matter. 

Results of the shear viscosity 7 and specific 
shear viscosity 7/s, where s is the entropy density, 
are shown in Fig.2. The temperature dependence of the 
shear viscosity is similar to that in Ref.[20] at different 
densities, while 7 increases with increasing density 
especially at lower temperatures due to the strong 
Pauli blocking effect. The specific shear viscosity 
decreases with increasing temperature, and it is similar 
in both magnitude and trend to those obtained from 
BUU calculations using the Green-Kubo formula’). It 
is interesting to see that at higher temperatures the 
specific shear viscosity is about 4~5 times the 


[11] 


lower limit from Ads/CFT calculation” , and this is 


already close to that of QGP extracted from the study 


12°! At lower 


using a viscous hydrodynamical mode 
temperatures the specific viscosity increases with 
increasing density due to the Pauli blocking effect, 
while at higher temperatures the dependence on the 
density is rather weak. 

Due to the sharper momentum distribution of 
neutrons compared to that of protons in neutron-rich 
nuclear matter, the total shear viscosity is dominated 
by neutrons which have a longer relaxation time in 
asymmetric nuclear matter compared to that in 
symmetric nuclear matter. This is confirmed in Fig.2 
that both the shear viscosity and the specific shear 
viscosity are larger in neutron-rich nuclear matter. In 
addition, it was seen!” that both the shear viscosity and 
satisfy the 


approximation with respect to the isospin asymmetry, 


specific shear viscosity parabolic 


i.e., 


m(p,T,6) = 9(p,T,6 =0)+ Mon (PT), (13) 


(2\io.r.0)=(2](o,7,0-04(2) (p,T)5?. (14) 


Similar to the symmetry energy, the second-order 
coefficient can thus be defined as the symmetry shear 
viscosity or the symmetry specific shear viscosity. The 
density and temperature dependence of them are 
shown in Fig.3. It is seen that both the symmetry shear 
viscosity and symmetry specific shear viscosity 
decrease with increasing temperature. At lower 
temperatures, both of them increase with increasing 
density. At higher 


dependence is rather weak. 7,,, and (1/5) ym are 


temperatures, the density 


important quantities in understanding transport 


properties of neutron-rich nuclear matter, and they 
deserve further studies in the future. 


era 300 E T T ta) (b) 300 
E \p =2 Pp L J 
> 200 \ 4 200 
£ à, MDI Interaction 
a 100 Fp Th sd + 100 
> T= 50 MeV ] 

+ a — 0 
z T=10MeV/ (d) g0 
E 4 60 
= ] 

g J9 
a 4 20 
Š r lE MOV] 

05 10 15 20 25 
T/MeV P-P" 


Fig.3 (Color online) Density and temperature dependence of 
symmetry shear viscosity ((a) and (b)) and symmetry specific 
shear viscosity ((c) and (d)). 


4 Conclusion 


Using a relaxation time approach, I studied the shear 


viscosity and specific shear viscosity of hot 


neutron-rich nuclear matter as that formed in 
intermediate-energy heavy-ion collisions by using an 
isospin- and momentum-dependent interaction. It is 
found that the specific shear viscosity decreases with 
with 
increasing density at lower temperatures due to the 
strong Pauli blocking effect. Furthermore, both the 


shear viscosity and specific shear viscosity are found 


increasing temperature, and it increases 


to increase with increasing isospin asymmetry of 
nuclear matter and roughly satisfy the parabolic 
approximation. The second-order coefficient in the 


050514-4 


XU Jun / Nuclear Science and Techniques 24 (2013) 050514 


expansion of the isospin asymmetry, which is defined 


as the symmetry shear viscosity or the symmetry 


specific shear viscosity, has also been studied. 
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